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' Abstract. We point out that, contrary to signs of heat capacities, thermodynamic 

Q . fluctuations are simply and unequivocally related to onset of instabilities that show up 

I near critical points. Fluctuation theory is then applied to Schwarzschild black holes 

surrounded by radiation. This shows that slowly evolving black holes along quasi- 
equilibrium states in cavities greater than 10^ Planck length will not evaporate below 
^ , the critical Hawking limit temperature despite the fact that pure radiation has a much 

00 ' higher entropy. 

m ' 
o 

In a mean field approximation, self-gravitating thermodynamic ensembles have local maxima of en- 
■rj" ', tropy or other Massieu functions derived from the entropy by a Legendre transformation, even though 
the potential energy is not bounded from below (see e. g. Horwitz and Katz 1978). These ensembles 



are generally not equivalent. This means that the domains of stability of different ensembles related 
I , by a Legendre transformation are different. Arguments based on signs of heat capacities and similar 
^j-jl thermodynamic coefficients are dubious criteria of stability. Far more secure results are derived from 
linear series {i. e. equilibrium curves) using Poincare's theorem about bifurcations or turning points 
(Ledoux 1958). One of the most elegant ways of using Poincare's method consists in following equi- 
^ I librium curves of pairs of conjugate parameters: inverse temperature versus energy, angular velocity 
d • versus angular momentum, and so on. (Katz 1978, 1979). This is the method used in Kaburaki et 
al. (1993) for Kerr black holes and in Katz et al. (1993) for Kerr-Newman black holes to find stabil- 
ity limits of stable configurations and degrees of instability in unstable configurations [see Parentani 
(1994) for the origin of the inequivalence and the degree of stability in canonical ensembles]. 

Here we want to stress that fluctuation theory is closely related to linear series of conjugate 
parameters and is indeed applicable to thermodynamic ensembles whether they are equivalent or 
nolj^ This is readily seen as follows. 

Think, for definition, of a Schwarzschild black hole in a sphere filled with photons at the same 
temperature (Hawking 1976). The total entropy of the system is S = S{E,V), where E is the total 
energy of hole and radiation and V is the volume of the cavity. Conjugate pairs with respect to S are 
(/3, E) and (P, V) with 

/3 = ^ P=l^ (1) 

^ dE' (3dV 
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Here f3~^ is the temperature, P the pressure (units are c = fi = G = k = l). Linear series of conjugate 
variables arc the equihbriiim curves I3{E) at fixed V (see figure 1) and P{V) at fixed E (not drawn). 
Changes of stabihty in the micro canonical ensemble (i. e. at fixed E and V) occur only at the critical 
point C, where there is a vertical tangent {{dp/dE)v = ±00). Positive slopes near C are on a branch 
of stable configurations, while negative slopes near C on one of unstable configurations. The proof 
is general and has been given elsewhere but is quite useful to remember and short enough to be 
reproduced here for the cases in which there is only one fiuctuating variable x. This is the situation 
of the Hawking example where x is the partition energy £'bh, the energy of the black hole and where 
the left over energy of the radiation is -Erad = E — -Ebh- 

Let S{x; E, V) be the entropy away from equilibrium. Equilibrium is defined by extremising the 
entropy with respect to x 

dj = (2) 

giving two solutions x = X{E, V) and x = Xi{E, V). Denote {X, Xi) collectively by Xa (a =nothing 
or 1) and S{Xa; E, V) = Sa{E, V). It is useful to define a temperature away from equilibrium 

^ = dES = Pix;E,V). (3) 

Clearly this function evaluated along the equilibrium configurations x = Xa gives 



(/3W.=/3„=(f|j^ (4) 



which is the usual equilibrium inverse temperature. The equilibrium curves (3{E) and Pi{E) are 
represented in figure 1. 

Stable equilibria exist if and only ii x = Xa is a configuration at which S is maximum: 

(dls)^ < 0. (5) 

Consider now equation (2) in which x is replaced by Xa; (2) reduces then to = and the following 
identity holds 



d{dj)a 

dE 



= {d,P)a + {dlS)a(^)^^0. (6) 



X 

V 



dE 

On the other hand, a further derivative of (4) with respect to E with V kept fixed gives 

fdPa\ (r,2 o\ . f» 3x fdXa 



Eliminating dXa/dE between (6) and (7) gives then 

Since changes of stability occur when {d'^S)a goes through zero, they manifest themselves by 
vertical slopes because for {d^S)a^O 



(9) 
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(see figure 1). The point of a vertical slope with coordinates (Ec, (3^) is a critical point. Near a critical 
point the slope is positive on the stable branch /9(-E) where {d^S)a < and negative on the unstable 
branch /3i (E) . 

Consider now S{x; E, V) near a point of stable equilibrium, at x = X + Ax: 

S = S + ^{dlS){Axf + O [{Axf] (10) 
or, according to (9) 

(Ag)^ (11) 

where Ap = ^ — (3 and 



> 0. (12) 

V 



Ap is a deviation ^{x) from stable equilibrium values P{E) induced by the fluctuation of x away 
from the equilibrium configuration X{E,V). The probability dW of such a fluctuation of f3{x) in a 
range (/3 + A(5,(3 + A/3 + d/3) is proportional to exp (5 — S) [see Landau and Lifshitz (1980)]. The 
normalization factor of dW is easy to find in this quadratic approximation: 



dW = = exp 

V27r(A/3)2 



l(A/3) 



2 (A/3)2 



dp (13) 



where 



(A/3)2= . (14) 



dE J 



V 



The mean square fluctuation (A/^)^ of is thus given near the critical point by the positive slope of 
f3{E,V), i. e. where the specific heat is negative. [Note however that in the canonical ensemble, at 
fixed /9, the mean square fluctuations of the total energy are given by the slope of —E{j3, V) as usual 
(see (19.6) in Callen 1985).] The real physical mean square fluctuations of temperatures in the bath 
of photons (A^rad)^ ^^nd in the black hole (A^bh)^ are neither equal to (A^)^ nor equal to each other 
because subsystems have different heat capacities. It is, however, not hard to show that near the 
critical point 



lim mSl! = lim = L (15) 

x^Xa (A^)2 ^^^a (A^)2 



Thus (A/3)2 as given in (14) or (9) represents, to a good degree of approximation, the quadratic 
fluctuations of temperature that diverge at a critical point. Similar fluctuations of a pressure function 
P{x;E,V) = 13-^dvS might have been considered. 

The more general relation between positive slopes and mean square fluctuations of a conjugate 
variable in linear series near the critical points where a change of stability takes place is obvious: if 
C (like P) is the conjugate function of Z which, like E, is a control parameter with respect to some 
Massieu function H(C,\ Z) (similar to S{x\ E, V)), the equilibrium value of ( is 
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We assume that second derivatives of H are bounded and that, at the critical point, {dxH)a 7^ 0. If 
{dxH)a = 0, vertical slopes may be transformed into bifurcations. But very small perturbation of H 
like H + exE, \e\ <C \xE/H\ will easily transform bifurcation into turning points (Thompson and Hut 
1973) with almost no change in stability limits. 

In these circumstances, near the critical point, the probability of ( fluctuating by = ( — is 
given by 



dW 



exp 



2(AC)2 



(AC)^ = ^>0. 



(17) 



When there are more variables x, we imagine that the x*'s have been taken in such a way that 
— {dfjH)a is diagonal, say SijXj. The \j are called the Poincare coefficients of stability. If the spectrum 
of Aj's is non-degenerate (which is always possible to achieve with a small perturbation of H), there 
will be a vertical slope at each point where a change of stability occurs and at the particular point 
d(^/dZ will be the same approximate value given by (9). 

For a simple application of fluctuation theory, we go back to (13) in the case of the Schwarzschild 
black hole in a cavity for which P{E), given in Hawking (1976) and Gibbons and Perry (1978), is 



E 



3 



(18) 



where L is the radius in Planck length unit. Suppose the black hole evolves through a series of quasi- 
equilibrium configurations by extracting slowly energy out of the box whose size is kept fixed. We 

know that at the critical point the hole will certainly evaporate. What we want to estimate is the 
probability per unit time that the black hole evapolates before reaching its stability limit as a result 
of thermal fluctuations. 
Prom (14) and (18) 



dEJv 



Stt 



/3c 



(19) 



in which 



1287r'' 
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(20) 



The probability of complete evaporation of the black hole is essentially the probability of the system 
fluctuating from P{E) to [3i{E). Indeed, as soon as [3 is slightly bigger than the probability of 
going back to j3 is negligible for 6'rad — Si ^ S — Si. Therefore, to estimate dW let us take for a 
typical value of A/3 the following one A/3 = f3 — (5^ (taking P — (5i would give an underestimate of 
that probability because the quadratic approximation is no longer valid over such an interval). With 
|(/3 — Pc)/Pc\ ^ 1 ^-^d with (18) and (19), dW as given in (13) becomes 



dW 




exp 




dp. 



(21) 



Hence, for a given value of Af3/f3^, the probability of evaporation depends only on the size of the box. 
For instance, say, A(5//3^ = 0.1 and L = 10^, dW = lO'^dp but for L = W^, dW = 8 ■ W'^^dp, i. e. 
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dW drops drastically for L > lO^. If A^//?^ = 0.01 and L = 10^, dW = 6 ■ but for L = 10^, 

dW c^2- 10-^°5d/3. 

Having found dW, we may now estimate the rate of evaporation. When the equilibrium configura- 
tion evolves slowly from P to C (see figure 1), the characteristic time scale of the return to equilibrium 
after a jump in energy of Ax is of the order of At ~ /3^Ax (see Zurek 1980). Since the characteristic 
time tc for a jump from point P at P{E) to point Q at Pi{E) is also proportional to it follows that 
the probability of evapolation per unit time is of the order of dW/ (3^(13^ — 13) (see Piran and Wald 
1982 for a similar argument). 

What the slope of [3{E) shows therefore is that for L > 10^, fluctuations of temperature leading to 
complete evapolation are totally negligible and the evolution of the black hole proceeds down almost 
to point C before evaporation takes place. The bigger the box, the closer to C. This holds in spite of 
the fact that near C, configurations are metastable and evaporated black holes (i. e. pure radiation) 
have considerably more entropy: 5rad — Sbh+rad ~ > 10^. 

Let us add that for L > 10^, backreaction and quantum gravity are totally negligible effects (York 
1985) and that the radiation is far from being general relativistic near the critical point (negligible 
self-gravity — see Parentani et al. 1994). Thus, Hawking's approximation for S on which our approxi- 
mations are based is perfectly good. Instead, for smaller boxes, the thermodynamic analysis probably 
loses its validity because the quantum fluctuations become important since one approaches the Planck 
temperature at point C. 
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Figure caption 

Figure 1. Qualitative features of Hawking's (3{E), Pi{E) and of f3j.g^{E). A/5(P) is the fluctuation 
around a point of equilibrium P. The critical fluctuation which will cause the evaporation of black hole 
arc of the order of (3{P) — l3i{Q). The vertical line E = Eb^V) is that for which S = S^adiE). For 
E < Eb the entropy S{E) is smaller than S,-i^{E) at fixed volume V. Thus, any point P between C 
and B represents a superheated metastable (in fact a stable -see the text) black hole in equilibrium 
with radiation. Correspondingly, the portion of E > Eb on the curve of l3j.g^{E) expresses Gibbons 
and Perry's (1978) superheated radiation state. 
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